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The problem of integer partitions is addressed using the microcanonical approach which is based on the analogy 
between this problem in the number theory and the calculation of microstates of a many-boson system. For 
ordinary (one-dimensional) partitions, the correction to the leading asymptotic is obtained. The estimate for the 
number of two-dimensional (plane) partitions coincides with known asymptotic results. 
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1. Introduction 

In this paper, we address the problem of integer partitions using the physical approach based on 
microcanonical treatment. Partitioning of integers is a problem in the number theory that originated 
in the works by Leibniz |1] and Euler [;2|. A partition of a positive integer n is a way of writing n as a 
sum of positive integers, where the order of the summands is insignificant. The number of partitions 
p(n) is called a partition function 1 3, Ch. 1] Further in this work we refer to p[n) simply as the number of 
partitions to avoid confusion with the respective physical term. 

To clarify the notion of partitions, let us consider the number 5. It can be represented as the following 
sums: 

5 = 4+1 = 3 + 2 = 3+1 + 1 = 2 + 2+1 = 2+1 + 1 + 1 = 1 + 1 + 1 + 1 + 1. 

Therefore, one has for the number of partitions p(5) = 7. 

Generalization for higher-dimensional partitions is made in the following fashion. In the D-dimensional 
case, an integer n is represented as a sum of positive integers _j B : 

n = E n h...i D > CU) 
h,...,i D »0 

where nt lm i D 3= rij 1 „j D whenever z'i s= j u i 2 =S j 2 , i D s= j D p. 179]. 

For instance, in the case of two-dimensional (plane) partitions of 3, one has the following possibilities: 

3 2+1 2 1+1+1 1+1 1 

+ + + 

1 1 1 (1.2) 

+ 
1 

yielding the number of plane partitions p2(3) = 6 . This can be easily shown by considering the 
numbers ; 2 to be elements of square lxl, 2x2,3x3, and so on, matrices. To obtain the sum of matrix 
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elements equal to 3 (provided that they are sorted in a non-ascending order in both rows and columns, 
according to the above definition), we have: 

/ 1 1 1 \ /10(M 

»(-)■(;:)■ CiM: ;;)■(:::)• 

Note that zero elements are skipped when writing multidimensional partitions. 

2. Physical analogy 

There is a straight analogy between the number of partitions and the number of microstates of a 
many-boson system of harmonic oscillators. The approach based on such an analogy was suggested as 
early as in 19S0 by Nanda |S], whose paper is now frequently overlooked, and was later utilized by a 
number of authors to solve related problems 

Hi. 

Let the system of one-dimensional harmonic oscillators with spectrum £■• = hw j have an energy E: 

E = ha)^jt, (2.1) 

i 

where jj denotes the quantum number of the zth particle. The summation runs only over the excited 
states and thus the particles in the ground state (with zero energy) can be arbitrary in number. The 
set {71 , 72. • • •} corresponds to a certain microstate of the system. In the quantum case, the particles are 
indistinguishable, so the permutation of {71, J2. ■ - ■} does not lead to a new microstate. At this point, we 
can match the set {71,72, •••} and a partition of the number n = E/hco. Thus, the number of microstates 
T(E) is equal to the number of partitions of n. In the set {71,72, . • ■}, any of the numbers jj can be equal, so 
one must consider bosonic oscillators. 

To make an asymptotic estimation of the number of microstates T(E), one can use the well-known 
Hardy-Ramanujan formula |9] for integer partitions: 

p HR ( „ ) = _L *v^73v^ (22) 
4V3« 

therefore, 

T(E)= " e nV2l3VEIf^_ (23) 

A\f?>Elh(D 

In the next section, we will derive this expression from physical considerations and obtain the first cor- 
rection to this asymptotic result. For simplicity, we further put the unit of energy hw = 1. 

The above considerations might be extended - with minor reservations - to higher-dimensional par- 
titions, the two-dimensional (plane) ones being best studied. The respective analysis is presented in sec- 
tiong] 

3. Corrections to the leading asymptotics of p{ri) 

Partition function Z{p) can be expressed as an integral by introducing the number of microstates 
T{E): 

OO 

Z{/3) = Y,e~ pEl '= f T{E)e- pE dE. 
I 

The above expression is nothing but the Laplace transform, so, by inverting it, we obtain: 

7+100 

r(JB)=-^T f ZfflePdp. (3.1) 
2m J 

y-ioo 
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The entropy S{/3) is equal to 



Thus, we have 



and 



SQ6) = fiE + In ZQ6). 



S(P) P -PE 



7+1CXD 



T{E) = — ( e m W. 
2m J 



y-loo 



Using the Taylor series for entropy in the vicinity of Pq: 



S{p) - S{Po) + j } S" (j8 ) (P-Po) 2 + ^ S'" (Po)(P-pof, 



where Pq is the stationary point, 



for the number of microstates, one obtains: 

y+ioo 

exp 



res) 



3 S(/3 ) 



2m 



s 



S'{p ) = 0, 



^s"(p )(p-p f + ^s'"ffioW - Pof 



dp. 



7-100 



Using the replacement p = ix+ /3q, we get: 

e S(/3 ) 



TIE) =* 



2n 

27T 

3 SC/3 ) 
27T 



f exp[-^S"(A)^ 2 -^S"'(A)ix : 

-CO 
CO 

J e -^S"(p )x* \ ' ,.3 

-OO 
CO 

J e -^S"[/3 )x 2 



dx 



S^o) 3 
3! X 



d.v 



S'"(/3 ) 3 
— x 



This integral can be expressed via the modified Bessel function of the second kind, namely: 

e^°> 2S"(p ) \ [S"(/3 )] 3 I ( [S"(A))] 3 ^ 

1 J " 2^ v^|S"'W| eXP t3[S"'W] 2 J 1/3 UtS'"^)] 2 ) - 

To calculate the entropy S(j6), we write the partition function of the oscillator system as follows: 

Z{p)=ll[l-e-Pi) , lnZ()S) = -£ln 1-e-W). 
, ; 1 ' ./' 1 

Using the Euler-Maclaurin formula to calculate the sum, one obtains 0]: 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



n 2 1 



1 



S{p) = PE + \aZ{p) = PE + — + - ln/3 - - ln(2jr) + . 

bp I 2 

Limiting ourselves to the first two terms, the stationary point yields: 

71 



Po = 



(3.9) 



(3.10) 



(3.11) 



(3.12) 
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With the same accuracy, we have: 



S"QS )=— E 312 , s"\Po)=~E 2 

71 



Collecting the whole expression (3.9) together, for the number of microstates we obtain 



i 



18^6£ 3/4 



exp 



28 2E 
27 V ~3~ 



Km 



1 Ize" 



II 71 V 3 



(3.13) 



(3.14) 



From this formula, there follows a correction to the main asymptotics of the number of partitions, as one 
substitutes E with n: 



p(n) ■ 



■ exp 



28 2n 
27 V T 



Km 



1 2n 
27 vT 



18\/6n 3/4 

Taking into account the asymptotic series expansion for large arguments jl 



/7i e 
2 sfz 



1 + I - 

z 



we immediately arrive at the leading asymptotic behavior given by equation (2T2 



(3.15) 



(3.16) 



4. Asymptotic behavior of plane partitions 

We consider the energy spectrum of a 2-dimensional system in the following form: 

e0'i.j2)=;i+j2. 



(4.1) 



To obtain T{E), we repeat the derivation presented in the previous section following |7], see also Hill . 
Thus, the entropy S{p) is equal to 



SQ6) = /3E + In ZQ6), 
where for energy spectrum (4.1) , the partition function is 

z(0= n (i-e-^) . n n [i-e" Wi+j2) ] . 

all the energies /i=l 7*2=1 

For the logarithm of the partition function Z(/3), one has: 

oo oo r . oo r , 

In Z(0 = - L L In [ 1 - ^ ] = - £ gy In [l - e~« ] , 

71 = 1 J2 = l ' j = l 

where the jth level degeneracy is equal to 



therefore, 



oo 

lnZ(j3) = -£;ln 1-e 
;=i 



•A/ 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



After applying the Euler-Maclaurin summation formula, the entropy can be expressed in such a form: 



((3) 1 



1 



S08) = PE + lnZ(/3) = PE + +— In p -- . 



(4.7) 
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The standard saddle-point method, for T{E) from equation H3.lt . yields 



exp[SQ3 )] 

T{E) = F H . (4.8) 
y/2nS"[p ) 



As in the previous section, for the stationary point /3q, so that S'ifio) = 0, one obtains 



and 



2<T(3) ll/3 



3 



(4.9) 



S"(A)) = ; - 1/ , £ 4/3 - (4.10) 
[2^(3)1 1/3 



Thus, the number of microstates is 



T(E) = — L= [2<-(3)] 7/36 £" 25/36 exp j - [2< (3)] 1/3 £ 2/3 - -\ , (4.11) 



67T 

Substituting energy E with integer we immediately obtain the result for plane partitions in the follow- 
ing way: 

p 2 in) = — L [2<T(3)] 7/36 n~ 25136 exp J § [2f(3)] 1/3 « 2/3 - H . (4.12) 
V^tt 12 6) 

Our estimation differs from the result of Wright 1 12] 



p^(n) = [2<T(3)] 7/36 n~ 25136 exp 1 - [2f (3)] 1/3 « 2/3 + c }• (4.13) 

VD7T I 2 



by a constant factor: 



c = t(-l) = -0.165421... versus - i = -0.166666... 

6 



Note that Wright's asymptotics was later confirmed by Nanda [fj with a much more sophisticated analysis 
of the oscillator system as compared to our approach. 



5. Results and discussion 

We performed calculations for the number of ordinary partitions pin) and plane partitions p^in) 
using the expressions d3.15> and <4.12) obtained in this work. 

The comparison with exact values and the leading asymptotics by Hardy and Ramanujan are given in 
tableland in ngureQ] In the figure, the relative errors are plotted. For n — 32 and n > 33, our correction 
becomes negative. For n > 20, the error is less than one per cent, and it never exceeds seven per cent 
except for n = l, where pin) = 1, and our formula yields 2, and n = 5 with our value of 8 versus pin) = 5, 
which leads to about 15 per cent error. Beyond this domain, the maximum relative deviation of 0.0083285 
is found for n = 186. 

The plane partitions are demonstrated in tableland in figure [2] Our result for plane partitions pro- 
vides a better estimation compared to that of Wright up to n = 7573. At n = 2679, expression <4.12> starts 
to underestimate the number of plane partitions while <4.13> asymptotically approaches the real values 
of p^iri) from above for all n. Interestingly, the approach to the calculation of plain partitions used in 
this work which takes into account the level degeneracy provides a much better estimate compared to an 
earlier suggested | IS] direct treatment of multidimensional oscillators which failed to produce a correct 
pre-exponential behavior of the power of n. 

To summarize, the expressions for the number of integer partitions were obtained from the analogy 
between this number-theoretical problem and the physical problem of calculating the states in a many- 
boson system within the microcanonical approach. The correction to the main asymptotics for ordinary 
(one-dimensional) partitions is shown to give a good estimate even for small numbers. 
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Table 1. Number of integer partitions. Real data p{n) are calculated by wxMaxima 0.8.7. Our result 13.15) is denoted as p om {n). 



TL 


pin) 




P (it) 


l 




1 


2 


2 


2 




? 


3 


2 


3 




Q 
O 


4 


3 


4 




c. 

J 


6 


5 


5 




7 


9 


8 


6 




1 1 


13 


11 


7 




1 R 


18 


16 


8 




99 


26 


23 


9 




ou 


35 


31 


10 




42 


48 


43 


20 






692 


631 


30 






6080 


5607 


40 




37338 


40080 


37244 


50 




204226 


217590 


203334 


60 




966467 


1024004 


961084 


70 




4087968 


4312670 


4061899 


80 




15796476 


16606782 


15686810 


90 




56634173 


59367760 


56218131 


100 




190569292 


199280893 


189113660 


150 




40853235313 


42369336269 


40515857434 


200 




3972999029388 


4100251432188 


3939941762556 


300 




9253082936723602 


9494094811675004 


9178612996544433 


400 




6727090051741041926 


6878471626940064454 


6675403544180355365 


500 


2 300 165 032 574 323 995 027 


2346386625611060168255 


2 283 287 889 193 750 426 745 


600 


458 004 788 008 144 308 553 622 


466396419561000383265349 


454 786 235 887 79 1 524 386 094 


/UU 


60 378285202834 474 611 028659 


Cl Af\1 CO/I 1 OC IOC (1(1(1 007 OCC OOO 

ol 4U1 334 l3bZi50U99 H5 1 000(592 


59 9 lU 590 l5Z 043 1(59 <L6a 1 5Z 2o5 


800 


5 733 052 172321 422 504456 91 1 979 


5 823 869 045 997 298 219 672 106 135 


5 695 741 93 1 526 069 206 199 862 70 1 


900 


415 873 68 1 190 459 054 784 1 14 365 430 


422 080 9 1 1 932 43 1 823 414 68 1 187 746 


413 258 060 326 350 156 090 055 206 82 1 


1000 


24 061 467 864 032 622 473 692 149 727 99 1 


24401996316802476288263414943062 


23 9 14 862 944 527 589 687 173 572 605 174 


10000 


3.61673 x 10 106 




3.63281 x 10 106 


3.60761 x 10 106 


100000 


2.74935 x 10 346 




2.75321 x 10 346 


2.74695 x 10 346 


200000 


1.14215 x 10 492 




1.14328 x 10 492 


1.14144 x 10 492 


500000 


1.52473 x 10 781 




1.52568 x 10 781 


1.52412 x 10 781 


1000000 


1.47168 x 10 1107 




1.47234 x 10 1107 


1.47127 x 10 1107 



Table 2. Number of plane integer partitions, real data pzin) are presented according to (l^Q. Our result is denoted as p° m {n). 



n 




pf(ri) 


p° m {n) 


1 


1 


2 


2 


2 


3 


3 


3 


3 


6 


7 


7 


4 


13 


14 


14 


5 


24 


27 


27 


6 


48 


51 


51 


7 


86 


94 


93 


8 


160 


169 


169 


9 


282 


300 


300 


10 


500 


526 


525 


15 


6879 


7174 


7165 


20 


75278 


77828 


77731 


25 


696033 


716466 


715574 


ou 


O OOO jDJ 




O OU / \JZJ 1 


40 


281846923 


287805195 


287446950 


50 


10499640707 


10690144561 


10676838029 


60 


314689799781 


319730938404 


319332954078 


70 


7937771067795 


8052251247267 


8042228226656 


80 


173781688194937 


176070453255058 


175851289976638 


90 


3376508618954817 


3417560875691699 


3413306874910366 


100 


59206066030052023 


59876276156314650 


59801745303809873 


125 


52658376905566496345 


53170908672841062211 


53104724310582212659 


150 


30 669 139 297 980 503 425 545 


30933055614865435615600 


30 894 551 767 372 734 822 993 


175 


12797410813969092203 145839 


12 896 659 9 14 533 283 896 907 383 


12 880 606 827 741 697 065 333 1 90 


200 


4066263490068623016919082185 


4 095 085 410 386 804 83 1 530 5 14 486 


4 089 988 062 550 422 858 47 1 57 1 83 1 


1000 


3.54259 x 10 84 


3.55112 x 10 84 


3.54670 x 10 84 


10000 


4.50750 x 10 401 


4.50983 x 10 401 


4.50422 x 10 401 


100000 


1.11796 x 10 1876 


1.11808 x 10 1876 


1.11669 x 10 1876 
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Figure 1. Comparison of relative errors for the estimations of the number of integer partitions from the 
real values (calculated by wxMaxima 0.8.7). Solid curve - our result I3.1S) . dashed curve - the leading 
asymptotics provided by the Hardy-Ramanujan formula (2.2) . Large n domain in an enlarged view is 
presented on the right. 



0.02 



iff 



— 0.01 



3_ 




Figure 2. Comparison of relative errors for the estimations of the number of plane partitions from the 
real values p2 in) (l^Q- Solid curve - our result (4.12) , dashed curve - Wright's formula (4.13) . 
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AcuMnTOTMMHi <|>opMy/ii/i fl/ifl po36i/iTTiB Lji.ni/ix m/ice/i y niflxofli 
MiKpoKaHOHiHHoro aHcaivi6/iK) 

fl,. npoxopoB, A. PoBeH L iaK 

Kacfieflpa TeopeTUHHOi (t>i3WKH, ZlbBiBCbKnti HaujOHa/ibHufi yHiBepcmeT iMeHi iBaHa <t>pam<a, 
Byn. flparoMaHOBa, 12, /lbBiB, 79005, YKpaiHa 

P03r/WHyTO 3afla4y npo po36nTra u,\nvtx 4wce/i y Mexax MiKpoKaHOHHHOro niflxofly, akhw rpyHTyerbCH Ha aHa- 
zioriY Mix L^icK) 3afla4eio 3 TeopiT 4wce/i i o64wc/ieHHflM Ki/ibKOcri MiKpociaHiB 6araTo6o30HHOi cwcreMH. fl/in 

3BW4aWHWX (OflHOBWMipHMX) p036nTTiB OTpWMaHO nOnpaBKy flO rO/lOBHOi' aCMMnTOTHKH. Ol^iHKa Ki/lbKOCTi flBO- 

BMMipHnx (n/iocKwx) po36wTTiB flo6pe y3rofl>KyeTbC5i 3 BiflOMUMH acuMnTOTUHHUMn pe3y^bTaTaiviw. 
K/i H) h o b i c/ioBa: po36tmR ujnwx mAcen, n/iocKi po36vnrn, 6o3ohhi cviaeMiA 
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